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ABSTRACT 


A study is made of post-buckling response and first-ply failure of flat unstiffened fibrous 
composite laminates subjected to uni-axial compression and a uniform temperature field. 
The analysis is based on the first-order shear deformation von-Karman type nonlinear 
plate theory. An efficient finite element formulation based on the principle of minimum 
potential energy is used for determining response characteristics, in conjunction with the 
first-ply failure analysis based on Tsai-Wu criterion. The main task is to obtain the 
stability boundary and the first-ply failure envelope for symmetric laminates having five 
different types of lay-up. The influence of temperature dependent thermo-mechanical 
material properties on the response of laminates is also investigated. 



CONTENTS 


Acknowledgements j v 

Abstract v 

Contents ' v j 

List of Figures viii 

List of Tables x 

Nomenclature x j 

1 Introduction 1 

1.1 Composite Materials 1 

1 .2 Literature Review 3 

1.3 Aim of Present Investigation 5 

2 Finite Element Formulation and Failure Model 7 

2.1 Introduction 7 

2.2 Plate Deformation, Strains and Stress-Strain Relations 7 

2.3 Stress Resultants at the Mid-Plane of the Laminate 9 

2.4 Total Langrangian and Virtual Work Equation 11 

2.5 Finite Element Formulation 13 

2.5.1 First- variation of total in-plane strain 14 

2.5.2 First-variation of bending strain 16 

2.5.3 First-variation shear strain 17 

2.5.4 Equivalent nodal load vector 17 

2.5.5 Non-linear equilibrium equation 18 



2.6 Solution of Non-linear Equilibrium Equation 18 

2.7 Failure Criteria 20 

3 Numerical Results and Discussion 22 

3.1 Convergence and Validation Study 22 

i 

3.2 Analysis Details 23 

3.3 Numerical Studies 27 

3.3.1 Thermal post-buckling response 27 

3.3.2 Stability envelope 28 

3.3.3 First-ply failure envelope 29 

3.4 Parametric study: effect of number of layers 37 

4 Concluding Remarks 39 

References 41 

Appendix I 45 

Appendix II 47 



LIST OF FIGURES 


2. 1 Geometry of a NL layered laminate 9 

2.2 Edge traction 12 

3 . 1 Convergence and validation of solution 23 

3.2 Verification of solution under combined loading 24 

3.3 Details of geometry and boundary conditions and sign 

conventions for stress resultants and generalized displacement 25 

3.4 Element and node numbering scheme for laminated plates 27 

3.5 Response of (±45/0/90)2 3 laminate considering thermal 

dependent (TD) and independent (T1D) material properties 28 

3.6 Stability envelope for (±45/0/90)2., laminate considering 

thermal dependent (TD) and independent (TID) material properties 29 

3.7 Stability envelopes considering thermal dependent (TD) and 

independent (TID) material properties 30 

3.8 Stability and failure envelopes for (+45/-45/0/90)2 S laminate 
considering thermal dependent (TD) and independent (TID) 

material properties 30 

3.9 First-ply failure envelopes considering thermal 

dependent (TD) and independent (TID) material properties 31 

3.10 Transverse deflection for (45/-45/0/90)2 S under 

combined loading at different temperature interval 32 

3.1 la Edge displacement under axial compression at x = 0, node 56 33 



3.11b Edge displacement under temperature rise at x = 0, node 56 33 

3.12a Variation of stress cr x at node 1, layer 2 34 

3.12b Variation of stress Txy at node 1, layer 2 35 

3. 12c Variation of stresses under temperature rise at node 1, layer 2 35 

3.13a Variation of stress a x at node 1, layer 15 36 

3.13b Variation of stress Txy at node 1, layer 15 36 

3.13c Variation of stresses under temperature rise at node 1, layer 15 37 

3.14 First-ply failure envelopes for (45/-45/0/90) ns laminates 

under combined loading 38 


LIST OF TABLES 

3.1 Composite panels considered 24 

3.2 Material mechanical properties 26 

3.3 Material strength properties ' 26 



NOMENCLATURE 


a, b, h 

C,j 

E, 

G,j 

K t 

k l ,k nl ,k„ 

ki, A'2 

M x .MyMxy 

b/1 x , M v , bA xy 

N 

NL 

N„ (1 ■1,N) 

Nx.Ny.Nxy 

N*,Nx,Nv 

Pn. ft 
Px.py 

Qx, Qy 

<? 

R 

R' 


Dimensions of the plate along x, y and z coordinate directions 
Elastic constants of a lamina in the principal material directions 
Young's moduli along the principal material directions 
Shear moduli in three orthogonal principal planes 
Tangent stiffness matrix 

Linear, nonlinear and initial component matrices of tangent 

stiffness matrix 

Shear correction factors 

Bending stress resultants per unit length 

Thermal moment resultants per unit length 
Number of nodes per element 
Number of layers of laminated plate 
Shape function 

In-plane stress resultants per unit length 
Thermal force resultants per unit length 
Local edge tractions per unit area 
Global edge tractions per unit area 
Transverse shear stress resultants 
Transverse force per unit area 
Equivalent nodal load vector 
Shear strength of the lamina in plane 2-3 



Vl2, V]3, V23 


AT 

{c> 



W 

L°) 

W 



vip), 

AM 

M 

W 

M 

M 

*0 




Poisson's ratios in planes 1-2, 1-3 and 2-3 , respectively 

Rise in temperature 

Bending strain vector 

Displacement of the i A node 

In-plane strain vector 

Transverse shear strain vector 

Linear component of in-palne strain vector 

Nonlinear component of in-palne strain vector 

Residual (unbalanced) force vector 

Incremental displacements at the end of the i* iteration 

Extensional matrix 

Extension-bending coupling matrix 

Bending stiffness matrix 

Transverse shear stiffness matrix 

First-variation of a quantity 

Partial derivative of a variable with respect to x 
Partial derivative of a variable with respect to y 

Partial derivative of a variable with respect to z 



s 


T 


T 


U, V, w 


UO, Vo. Wo 


Xu x c 


Yu Yc 


Z(, Zr 


a 

a y 


cti,a 2 


0 X , Oy 


<^X, CTy, Txy 


CTj, OS, (73 


VXZ. Tyt 


?n, t] 3 , T23 


Shear strength of the lamina in plane 1-3 
Shear strength of the lamina in plane 1-2 
Temperature at the midplane of the plate 

Displacement along x, y, z respectively of any point (x, y, z) in the 
plate 

Mid-plane displacement in x, y and Z directions, respectively 
Tensile and compressive strength, respectively; of lamina in the 
fiber direction 

Tensile and compressive strength, respectively, of lamina in the 
direction transverse to the fiber direction 

Tensile and compressive strength, respectively, of lamina in 
principal material direction-3 

Fiber orientation measured +ve in anti-clockwise direction 
Coefficients of thermal expansion with respect to the structural 
axis 

Coefficients of thermal expansion of individual layers in direction 

of fibers and normal to fiber direction respectively 

Rotation of normal to the undeformed mid-plane in xz and yz 

planes, respectively 

Stress components in xy plane 

Normal stress components in principal material directions 1,2 and 
3, respectively 

Transverse shear stress components in xz and yz planes, 
respectively 

Shear stress components in planes 1-2, 1-3 and 2-3 respectively 


CHAPTER 1 


INTRODUCTION 


1.1 Composite Material 

Engineering sounds metallic. Metals are, by far, the most versatile engineering materials. 

In modem engineering structures, new composite material systems exhibiting exotic 

\ 

properties play a big role as compared to the traditional metallic ones. 

The word ’’composite" refers to a material consisting of two or more chemically distinct 
constituents, on a macro scale, having a distinct interface separating them. In addition, 
the properties are strongly influenced by the properties of its constituent materials, their 
distribution and their mutual interaction. This definition encompasses the 'continuous 
fiber-reinforced laminated composites' which are of primary interest in this study. 

A laminated composite comprises several layers, termed lamina, consisting of high 
performance fibers (e.g. glass, boron, graphite) unified by advanced binders called matrix 
(e.g. epoxy resin). The use of composite material in various structural applications is 
dictated by factors like the increased strength, stiffness, improved structural response, 
low maintenance cost and high service life. Because of their high strength and stiffness to 
weight ratios and better thermal properties, these new materials are going to play a great 
role in the construction of aeronautical and aerospace vehicles as well as of turbojet 
engines, large space structures and nuclear reactor components. 

Composite structures are generally subjected to external mechanical loads, aerodynamic 
loads, gust loads, thermal loads, moisture loads or its combinations. A major part of the 
above mentioned structures are going to operate in hostile thermo-mechanical 
environments which could result in significant modifications in their response 




characteristics. This is especially valid for supersonic space vehicles where, as a result of 
aerodynamic heating, the structure is exposed to a non-uniform temperature field. In 
view of this, the consideration of simultaneous actions of thermal and mechanical loads 
pose challenging technical problems in understanding the structural response, failure 
criteria and design of composite structures. 

High performance composite laminates lead to new opportunities in the design of 
structures. Composite structures show flexible anisotropy and are characterized by 
bending-extension coupling. This gives the freedom to tailor the stiffness and strength, 
especially in the case of laminated filamentary composites. In contrast to an isotropic 
material, composite materials have a low ratio of transverse shear modulus to the in- 
plane modulus. Owing to this fact, the transverse shear deformation plays a much more 
important role in reducing the effective flexural stiffness of laminated composites plates 
than in the corresponding metallic plates. The laminates used in most of the structures are 
usually thin and in many applications, the design may require an exhaustive use of the 
post-buckling strength. In the presence of elevated temperature the response 
characteristics become more complicated. Further, in many composite materials the 
temperature causes significant variations in the thermal and mechanical properties. The 
strength and moduli of most of the materials decrease with rise in temperature, while 
Poisson's ratios and coefficients of thermal expansions usually increase. 

Design has to make use of the above mentioned features to provide better performance 
advantages depending upon the functional requirements of the structure. This demands 
an efficient modeling. 

A lamina is idealized as an orthotropic material. Plates and shell structures made of 
laminated composite materials are often modeled as equivalent two dimensional 
continuum using classical laminate theory, which is a direct extension of the Kirchoff 
hypothesis for homogeneous thin plates. As is well known, geometric non-linearities 
arise from finite deformations of an elastic body. The investigations of the response in 
bending, buckling and post-buckling of composite laminates is far more involved and 
complicated when compared to its isotropic counterpart. The von-Karman's assumptions 
governing large transverse deflections (of the order of thickness) for anisotropic, 
heterogeneous, thermoelastic thin plates are extended for the non-linear analysis of 



laminated composites. For moderately thick plates, in which the effect of shear 
deformation can be significant, the application of the Reissner-Mindlin theory or higher 
order shear deformation theories becomes necessary. 

It is almost impossible to look for closed form solutions for most of the practical 
problems involving composite laminates. The approximate analytical techniques, which 
have been applied to study the buckling and the post-buckling behavior of flat composite 
plates, are mostly limited to the cases having simple boundary conditions. Using the 
finite element technique, it has been possible to perform the analysis of composite 
laminates of complex geometry (e.g. plates with cutouts or having stiffeners), with 
complicated boundary conditions (e.g. elastic, discontinuous or point supports), having 
variable thickness and with temperature dependent material properties. The displacement 
finite element models, based on the principles of virtual displacement, have emerged as 
the most powerful tools for the solution of problems related to composite structures. 

An efficient and economical design of laminates by fully exploiting their strength 
characteristics with high reliability and safety requires a thorough understanding of the 
failure mechanism. A first ply failure analysis is often used in predicting the laminate 
strength. Various macroscopic failure criteria based on tensile, compressive and shear 
strength of the individual lamina are used for this purpose. Equally important is to have 
the knowledge of the mode of failure and its consideration. Laminated composites may 
fail by fiber yielding, matrix yielding, fiber breakage, delamination of layer or by 
fracture. 

In case of thermal problems, the elevated temperature usually lowers the matrix strength 
such as transverse and shear strengths, but reduces residual stresses. The static tensile 
strength of laminates with fiber dominant failure modes is relatively insensitive to 
temperature, variations. 

1.2 Literature Review 

A vast amount of literature is available on diverse aspects of thermoelastic response of 
composite laminates. A comprehensive summary of the state of art in nonlinear analysis 



of bending and post-buckling of isotropic, orthotropic, anisotropic and unsymmetrically 
laminated plates was compiled by Chia (1988). 

A survey paper by Tauchert (1991) on thermally induced flexure, buckling and vibration 
of plates gives an up-to-date account of buckling, post-buckling and large deformation of 
plates associated with elevated temperatures. The most comprehensive review paper by 
Noor and Burton (1991) details hierarchy of composite models, the effect of temperature 
dependence of the material properties on the response, and the sensitivity of thermo- 
mechanical load to the variations in material parameters. 

During the last decade, many investigations have been done on the post buckling 
response of composite laminates. Huang and Tauchert (1988a, 1988b) investigated the 
thermal post-buckling of antisymmetric angle-ply laminated plates under uniform and 
parabolic temperature fields through a direct applications of the principle of minimum 
potential energy. A major effort in dealing with the problem of thermally induced large 
deflections by the finite element method was done by Chen and Chen (1989). Huang and 
Tauchert (1991) investigated the large deflection response of flat plates and cylindrical 
panels under thermal loading, using the finite element formulation based on the first 
order shear deformation theory. The nonlinear response of antisymmetric cross-ply 
laminates under the uniform temperature field was studied by Kumar and Dixit (1994). 

Chen and Chen (1991) examined the effect of temperature dependent properties on 
thermal post-buckling behavior of laminated plates. It was found that the post-buckling 
strength got reduced significantly when temperature dependent properties were taken into 
consideration. Mayers and Hayer (1992) studied the thermal post-buckling of 
symmetrically laminated plates using the Raleigh-Ritz procedure. Shen and Lin (1995) 
considered the post-buckling response of geometrically imperfect laminates using the 
Galerkin perturbation technique under uniform/nonuniform temperature loading. 

The buckling problem of isotropic long plates subjected to heat and axial loads was 
considered by Boley and Weiner (1960). Thermal buckling loads of initially compressed 
transversely isotropic and anti-symmetric cross-ply laminated thick plates were evaluated 
using the Galerkin method by Chen et al (1982). The influence of uniformly elevated 
temperature on the buckling loads of uniaxially compressed isotropic long plates was 
analyzed by Sadovsky (1993). Noor and Peters (1992) made an analytical study, based 



on first order shear deformation theory, for the thermo-mechanical buckling of composite 
plates subjected to combined thermal and axial loading. A mixed finite element method 
is used, for determining the stability boundary and for evaluating the sensitivity of the 
buckling response to variations in the material and geometric parameters. Noor et al. 
( 1993 ) calculated the buckling load and the post-buckling response of flat unstiffened 
panels subjected to combined temperature change and applied edge displacement. The 
analysis was based on a first order shear deformation, von-Karaman type nonlinear plate 
theory. Numerical results were presented showing the effects of variations in the 
laminate stacking sequence, fiber orientation, number of layers and aspect ratio of the 
panels on the thermo-mechanical buckling load and the post-buckling response. Liberscu 
and Souza (1993) analyzed the post-buckling of an imperfect, shear- deformable, 
transversely isotropic plate under combined thermal and compressive edge loading. A 
post-buckling analysis for a simply supported, composite laminated plate resting on a 
two-parameter elastic foundation subjected to combined axial compression and uniform 
temperature loading was conducted by Shen and Sun ( 1996 ). Argyris and Tenek ( 1995 ) 
studied the thermo-mechanical post-buckling response of fibrous composite laminates. 
Temperature-dependent material properties in cubic variations were incorporated. It was 
found that the mechanical property degradation due to temperature rise had large 
influence on the response of the laminate. Srikanth ( 1997 ) studied the response and the 
first-ply failure of thin temperature sensitive laminated composite plates under thermal 
loading. 

Quasi-static thermal stress response and failure behavior of graphite/epoxy laminated 
plate under combined thermal and mechanical loads was investigated by Chen et al. 
(1985). Effects of ablation, degradation of thermo-physical and strength properties at 
elevated temperature, and radiation and convective heat losses were included in the 
formulation. The maximum stress criterion was adopted for the failure predictions. 

1.3 Aim of Present Investigation 

The literature review reveals that very few investigations are available related to buckling 
and post-buckling behavior of composite plates subjected to combined thermal and 
mechanical loads. Also, studies related to the failure of composite plates under thermal 



loading are rather limited in scope. In much of the earlier works, an assumption is 
frequently made that the material properties are temperature independent. But, 
temperature causes significant variation in thermal, mechanical and strength properties of 
the material. The accurate prediction of response requires accounting for the effect of 
temperature in thermo-mechanical computational models. 

The present study focuses on having a better understanding of the buckling, post- 
buckling and failure characteristics of multilayered composite panels subjected to 
combined mechanical and thermal loads. Specifically, the objectives are 

• to obtain the stability boundary of different types of laminates under axial 
compression and uniform temperature loading. 

• to obtain the failure envelope of laminates in terms of the first-ply failure under 
combined thermal and axial loads. 

• to study the effect of temperature sensitivity of mechanical and strength properties on 
the response and failure of laminates. 

The study is performed using a finite element formulation based on the first-order shear 
deformation, von-Karman type non-linear plate theory. Most notably, temperature- 
dependent thermo-mechanical properties are also incorporated. 



CHAPTER 2 


FINITE ELEMENT FORMULATION AND 
FAILURE MODEL 


2.1 Introduction 

The aim of this chapter is to provide a geometrically nonlinear Mindlin type formulation 
for laminates. The finite element equations are derived based on the principle of 
minimum potential energy. The formulation is basically the one developed by Singh et al. 
(1996). It is based on Mindlin's plate theory (see, e.g., Hinton et al^ 1975) which accounts 
for the through-the-thickness shearing effects, resulting in, plane sections normal to the 
midsurface remaining plane but not normal to the midsurface after deformation. This 
model has the advantage of requiring only C° continuity for all the primary variables, and 
hence, enables an easy formulation. Such a model readily accommodates moderately 
thick plates; however, they can exhibit locking when applied to thin plates. Extensive 
studies on Mindlin plate elements (Paugh et al., 1978) envisage that the nine noded 
Lagrangian element, with selective integration of the shear terms, is quite adequate for 
predicting the response characteristics of both thin and thick plates. 

2.2 Plate Deformation, Strains and Stress-Strain Relations 

In the Mindlin Plate theory, the displacements u=[u, v, wj T at (x, y, z) are expressed as 
functions of the mid-plane ( z=0 ) translations u 0 , v 0 and m> 0 , and of the independent normal 
rotations &x and By in the xz and yz planes, respectively: 


u{x,y,z) = u 0 (x,y) + z$ x (x,y) 


( 2 . 1 ) 



v ( x >y,z) = v o (x,y) + z0 y (x,y) 
w(x, y,z) = w 0 (x,y) 

For a Mindlin plate, the relevant Green’s strain vector (Chia, 1980) is given as, 
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Introducing von Karman’s assumptions (Fung, 1965) and noting that w is independent of 


z, eqn. (2.4) can be rewritten, using eqns. (2. 1-2.3), as 
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where the suffix 7* and *nT stand for the linear and nonlinear parts, respectively, of the 
mid-plane strain components. 



Based on the assumption that a z is negligible, the thermo-mechanical constitutive 
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where k u k? arc the shear correction factors (Mindlin, 1951), Q tJ are the transformed 
reduced stiffness coefficients of a lamina (see Appendix I), a x ,a y ,a xy are the 

transformed thermal expansion coefficients and AT is the temperature change. 

2.3 Stress Resultants at the Mid-Plane of the Laminate 

The stress resultants for a laminate with £ NL’ number of layers as shown in fig.2.1 are 

given by 
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Fig. 2.1. Geometry of a NL layered laminate. 



Substituting eqn. (2.5) in eqn. (2.13) and integrating through the thickness, the expression 
for stress resultants become 
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2.4 Total Lagrangian and Virtual Work Equation 

The formulation is based on the virtual work equation for a continuum in the total 
Langrangian coordinate systems under assumptions of small strains and conservative 
loads. In the absence of body forces the virtual work equation can be written as, 

jSe T adJ'’ r -jSu r pdA = 0 (2.27) 

where V is the undeformed volume, a is the Piola-Kirchoff stress vector and 8s is 
variation of Green’s strain vector due to variations of displacement u (5 u) and p is the 
surface traction acting over an undeformed area A. Here 8 denotes the first-variation. The 

variation of strain energy is given by, 

8U = £ 8s r cnlV = £ + 8s y a y + 8y x> , + 8y V2 r y , + 8y a r xz ]dzdA (2.28) 


Substituting the strain expressions of eqns. (2. 5-2.8) and eqn. (2.1 1) in the eqn. (2.28) and 

using the stress resultant expression of eqns. (2. 14-2. 16), §U can be written as 

8U = f + 8c° y N y + Sy^N^ + 8K X M X + 8K y M y + SK^M^ + 8y xz Q x + 8y yx Q y ]dA 

(2.29) 

The externa! work done due to edge traction p„ and p t which are the local normal and the 
tangential tractions, respectively, acting on edge (see fig. 2.2) and written in terms of the 


global edge tractions p x and p y , is given by 
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The global edge traction _p*and p y are given in terms of local edge traction p n and pt as. 
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where <j> is the orientation of local normal traction p„ with respect to global direction x 
(fig. 2.2). Substituting for u and v from eqns. (2. 1-2.2) in the enq. (2.30) and carrying out 
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the integration through the thickness, the variation of external potential energy due to in- 
plane forces becomes 

<K = +Sv a P y )ds (2.32) 

where 

^j hi h 2 n (Px,P?h (2.33) 

The variation of external potential energy with surface force is given by 

= ~^8* 0 qdA (2.34) 

where q is the intensity of transverse loading. The variation of total potential energy, i.e., 
eqn. ( 2 . 27 ), can now be written in terms of the mid-plane forces and strain/displacement 

quantities as 
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2.5 Finite Element Formulation 


The displacement within an element is expressed as a function of discrete nodal 
displacement: 

W = k o x ,& y Y = Z [*, 3 k ) ( 2 - 36 ) 

where 

7V= number of nodes in an element, 

Ni, V /' e [ 1 , N ] are the shape functions. 

If is 5 X 5 unit matrix, 

[«* » v w , u' 0l , 0 „ , 9„ j 7 is the nodal displacement vector, the superscript V refers to 

an element and the superscript l T indicates transpose of a matrix. The shape functions of 
a nine-noded Lagrangian element are given in Zienkiewicz (1992). 


The virtual work equation (2.35) can be re-written as 
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k’ }“ xy £ 

{K)=\K„K y ,Kj 


(2.37) 
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(2.39) 

(2.40) 

(2.41) 


are, respectively, the in-plane, the bending and the transverse shear strain components of 
the total mid-plane strains of eqn. (2.7) and 




(2.42) 

(2.43) 



2-5.1 First - variation of total in-nlanc strain 


The in-plane strain vector ^ 0 p }~[^' x ,£ y ,y x> ,] T Is expressed as a sum of linear and 
nonlinear components as given below: 
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(2.46) 


Linear and nonlinear components of the in-plane strain can be written in terms of the 

elemental nodal displacement as, 

(2.47) 
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kY - jDUW 


(2.48) 


where 


1 = K.i . Bol.2 . J (2.49) 

] = lAw /,,1 > ft ON I.. 2 7 ftoNL,3 >•••> ft ONI, N J (2.50) 

Here, Bo/., and Bow., are the linear and the nonlinear strain-displacement matrices, 
respectively. The elemental nodal displacement vector is given as 

W s h. fl 2» a 3 ( 2 - 51 ) 


where 



is the nodal displacement vector. 


The component matrices of the elemental strain-displacement matrix are given as 
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The component matrices of G are 
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(2.57) 


The in-plane strain vector is now given in terms of Ae and the Cartesian derivatives of 

shape functions as 



(2.58) 


The first-variation of the total in-plane strain at the mid-surface is obtained by taking the 
first variation of eqn. (2.58): 

4 "}= fcPW GStf (2.59) 
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Expressing w in terms of [A 6 ] matrix as function of shape functions and re-arranging 
the first-variation of total in-plane strain can be expressed as 
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2.5.2 First-variation of bending strain 

The bending strain {K }= [fC x , K yt ] can be expressed 
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Substituting eqns. (2. 1-2.3) in eqn. (2.63), the bending strain is expressed in terms of 

shape function as 

KMtfJM' (2 64) 

where 

{a}* is the elemental nodal displacement vector as defined in eqn. (2.51). and 

(2.65) 

in which the component matrices are 
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( 2 . 66 ) 


The first-variation of the bending strain can now be expressed as 


(2.67) 



2.5.3 First-variation of shear strain 


The transverse shear strain, }= f 
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can be expressed in terms of shape function as 

where 
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in which the component matrices are 
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The first-variation of shear strain is now expressed in terms of shape function as 


(2.72) 


2.5.4 Equivalent nodal load vector 


The external virtual work due to in-plane loads and transverse forces in terms of shape 
function is expressed as 
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in which the component matrices are 



K/J=^[/ 3 ] (2.78) 

where N t is the shape function corresponding to the i A node and [/ 3 ] is a 5 X 5 unit 
matrix. 

2.5.5 Nonlinear equilibrium equation 

The total virtual work equation (2.35) is discretised for an element by substitution of 
eqns. (2.62), (2.67), (2.72) and (2.74); 

ft}+ ft] r £w}+ [aj &]d4 

-w T t k r ip. }<« + ( ft r ft }*]= o <Z79) 

from the equation. (2. 1 8) 

44T>»rMri«JM+ftrak 

-4>n>jft,MBJM]<« (2.80) 

- swr [(k r ft + £.[s„ r ft }*]= ° 

Since the nodal virtual displacement are arbitrary, the element nonlinear equilibrium 

equation becomes 

*i>h 1 Ik r &)* ft r Ph ft r m- tu r ft, U - * - o w 

where 

« = ft, ft ft r {M, U + 1 ft f ft }<« + (ft f ft ft ( 2 . 82 ) 

The eqn. (2.81) has a dual role of representing the total equilibrium equation; in an 
element form and in the assembled form, 

2.6 Solution of Nonlinear Equilibrium Equation 

The solution algorithm for the assembled nonlinear equations is based on the Newton- 
Raphson method. If an initial estimate of displacement vector gives the residual forces 
y/\i } * 0, then an improved solution {a} i+1 is obtained by equating to zero the linearized 

Taylor’s series expansion of in the neighbourhood of £*}, as 



^{*L =^H =0 (2.83) 

where Ay?} is the incremental displacement vector, and K r is known as the tangent 
stiffness matrix evaluated at £?} and is given by 


K 

r "L aW 


(2.84) 


The linearized approximation to the relation between residual forces and incremental 
displacements at a point in the equilibrium path is given by eqn (2.83). The improved 
solution is then found as 

(2.85) 


The eqns. (2.83) and (2.84) represent the Newton-Raphson technique to the solution of 
nonlinear eqn. (2.81), which requires repeated calculation, and inversion of K r matrix. 

To improve on the numerical stability and convergence of the solution, the load is applied 
in small increments. The iterative solution is checked for convergence using the 

following criteria. 

i 

x 1 00 <; /? (2.86) 

where /? is sufficiently small number. 
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The tangent stiffness matrix K r cm be re-written as 
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(2.87) 
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in which the subscript ‘m’ refers to the total degree of freedom per element. 

dR 

Differentiating yr with respect to a. and noting that -^- = 0 for conservative loads, the 
tangent stiffness matrix assumes the following form (see Appendix II). 
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Substituting for the matrices B 0 ,B b , B s and G in the preceding equation one obtains 

K T =K L +K m +K (r (2.89) 

where the constant linear stiffness matrix X). is 

a-,. -iWmiBcj+mmjfa 

JA n 90) 

+ J, k,I[*M]+ [BJM [b ( HbJ[£][b,]Jm 

the initial displacement matrix K^l is 
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and the initial stress matrix K„ is 


A'.-rw 


N 1 

T x v [G\dA 

N A 


(2.92) 


The integration of expression for '{'{a} and Kj- is carried out using the Gaussian 
quadrature. A selective integration scheme is adopted with 3X3 integration rule to 
evaluate integrals of the functions of the membrane and the bending behavior and a 2 X 2 
integration rule is used for the transverse shear component. Predominance of zeroes in 
matrices Bo, Bh, B s , and G indicates that the expression for ^{aj and K T are best 


calculated explicitly, giving only the non-zero values 


2.7 Failure Criteria 

A detailed review of various failure criteria has been made by Reddy&Pandey (1987) and 
Reddy & Reddy (1992). Herein, for predicting the laminate failure, the first-ply failure 
analysis is done based on the Tsai-Wu criterion (1971); it is postulated that a failure 
surface in stress space exists in the form 


FiOj + FijOiOj-l 


i,j = l,...,6 


(2.93) 



wherein Fj and Fj are strength tensors of the second and fourth rank, respectively, and the 
usual contracted stress notation is used except that a 4 = x 23 , cr 5 = x 3u and a 6 = x 12 (in the 

principal material direction). The tensor polynomial form of the Tsai-Wu criterion is 
given below : 

Fia, + F 2 ct 2 + F 3 CT 3 + 2F 12 cnCT 2 + 2F 13 cr 1 a 3 + 2F 23 a 2 o 3 + Fi,cvi 2 + F 22 a 2 2 + F 33 ct 3 2 
+ F44C4 2 + F55CJ5 2 + F 66 ct6 2 > 1 (2.94) 
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And all other strength tensor components are zero. 



For the composite laminates under plane stress conditions, lamina failure is said to have 
occurred when the state of stress at any point with in the lamina (at mid-thickness) 
satisfies the tensor polynomial form of the Tsai-Wu criterion in which terms associated 
with normal stress component in third principal material directions are omitted. Equation 

(2.94) becomes 

Fi<Ji + F 2 o 2 + 2Fi 2 <J)<t 2 + F„a, 2 + F 22 a 2 2 + F 44 04 2 + F 55 CT 5 2 + F 6606 2 S: 1 (2.95) 

The nodal point stresses are calculated using the finite element model and the stresses are 
checked for the above failure condition. The failure of a node in the lamina represents the 
failure of that particular layer. The first-ply failure load refers to the first instant at which 
one, or more than one ply fails at the same load. 



CHAPTER 3 


NUMERICAL RESULTS AND DISCUSSION 


In this chapter, numerical illustrations are given for symmetrically laminated plates 
subjected to combined thermal and compressive edge loading. The finite element 
formulation presented in the previous chapter is used to study the response, mainly in 
terms of the stability and the first-ply failure envelopes. Regarding the prediction of the 
first-ply failure, the Tsai-Wu failure criterion is incorporated in the computational model. 
Numerical studies deal with the cases of perfectly bonded symmetric Graphite/Epoxy 
laminates having temperature dependent/independent, mechanical properties. 

3.1 Convergence and Validation Study 

The capability of the finite element formulation in predicting the post-buckling response 
under thermal loads is demonstrated by comparing results with those obtained by Huang 
and Tauchert (1991). An infinite layered regular anti-symmetric angle-ply laminate, with 
temperature independent material properties and subjected to uniform temperature rise, is 
considered. The plate is square, simply supported with edges restrained against the 
inplane normal displacement and free to move in the tangential direction. Figure 3.1 
gives temperature versus central deflection curves as obtained by using successive 
refined meshes. It is observed that results of 5X5 and 6X6 meshes are in good agreement 
with those of Huang and Tauchert (1991). To save computing time without sacrificing 
accuracy, the 5X5 mesh is adopted in the entire analysis. In addition, the computer 
program has checked well with the results given by Chen and Chen (1989) and Noor et 
al (1993) for the thermal post-buckling response under different boundary conditions, 
laminate sequences and material properties. It may be mentioned here that the validation 




study for the stresses under inplane mechanical loads has been reported by Singh et al. 
(1997). 


Temperature, T (°C) x 10 



Wc^/h 
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(X T ■■9x10*. L, * L : =0.254m. b/h * 125, (9=30°) 

Fig. 3.1 Convergence and validation of 

solution 

The program is also checked for the case of combined thermal and mechanical loads by 
comparing the results of the buckling strength obtained by Noor and Peters (1992) for 
(+45/-45/0/90>j» and (+30/-30)4. laminates. The laminates are simply supported, with 
edges free to move in the normal direction but restrained in the tangential direction, and 
subjected to a compressive load and a uniform temperature rise. As can be seen from fig. 
3.2, the results match quite well. 

3.2 Analysis Details 

In the present study, five types of graphite/epoxy composites are considered; the 
lamination scheme is given in Table 3.1. It is assumed that the laminate has perfectly 




(hj — 130.3 GPa, I; i -9.37 GPa, Gj.r -4.502 GPa, Grr —1.724 GPa, vjj — 0.33, ai= 0.139xl0’ <s , a? 
-9xl0' ft , L, - I.j 0.254m. h/h --- 125) 


Hg. 3.2 Verification of solution under combined loading 

bonded layers and during loading no delamination occurs. The panels are subjected to a 
combined action of the axial compressive load, N x and the uniform temperature rise, T. 


Table 3.1 Composite panels considered 


lamination 

scheme 


(±45/0/90)2. 


(±45/0 2 )2. 


(±45/0 6 ). 


(±45) 4 . 


(0/90)4. 


The geometric properties and boundary conditions are given in fig. 3.3; these are the 
same as adopted by Moor and Peters (1992). The boundary conditions correspond to a 
simply supported laminate, with normal displacements permitted and shear 
displacements restrained on all edges. It may be noted that at comer points all the five 
degrees of freedom are restrained. The material properties of the laminates are taken as in 
Chen el at, (1985). Therein, temperature-dependent mechanical properties are derived 
analytically based on the rule of mixtures in conjunction with a glass transition 
temperature of 200°C for epoxy. In addition, the graphite fiber is assumed to retain its 
stiffness upto 3300°C. The elastic moduli and thermal expansion coefficients are listed in 











Table 3.2, in which E, G and /u represent the Young's modulus, the shear modulus and 
the Poisson s ratio, while ct denotes the coefficient of thermal expansion. 


Geometric parameters 

Lj ~ 254 mm 
L ; = 254 mm 
Thickness of individual 
layers=0.127 mm 


Boundary conditions 

At x = 0, a: 

v = w = 6 y = 0 

At y = 0, b 

u=w = = 0 



u 


Fig. 3.J Details of geometry, boundary conditions and sign conventions for stress resultants 

and generalized displacements. 


Regarding the strength, for temperature exceeding 200°C, it is assumed that Gr/Ep 
strength decreases linearly to zero at a temperature when epoxy has been completely 



cooked out. The net effect of the epoxy cooked out is to deplete the ability for the load 
transfer from one fiber strand to another. Based on the limited number of experiments on 
Gr/F.p at high temperature, the temperature-dependent strength properties of a single ply 
Gr/Ep composite were developed, and are tabulated in Table 3.3. (Chen et al 1985). 
Mechanical and the strength properties are assumed to have a linear variation for 
intermediate temperatures. In the analysis, the initial body temperature of the plate is 
taken as 20°C. It implies that at this temperature the laminate is stress free. 


Table 3.2: Material mechanical properties 
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Tabic 3.3: Material strength properties 
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As mentioned earlier, numerical results are obtained by using a 5X5 discretisation for the 
full laminate. The details of the finite element mesh and locations of Gauss points with in 
the elements are shown in fig. 3.4. Although results are obtained for all five lamination 





































schemes given in Table 3.1 ? the detailed investigation, whenever required, is carried out 
for the quasi-isotropic laminate only. 
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Fig. 3.4 Element and node numbering scheme for the laminated plate 


3,3 Numerical Studies 

3.3.1 Thermal post-buckling response 

The thermal post-buckling response of 16-layered quasi-isotropic panels under uniform 
temperature loading is shown in fig. 3.5. There is a remarkable variation in the response 
and failure temperature for thermal dependent and independent material properties. In the 
process of temperature rise, temperature dependent properties reduce the stiffness of the 
laminate, thus causing variations in pre and post buckling deformations, critical 
temperature etc, 

For Graphite T- poxy composites, the response reveals a snap-through type behavior on 
the application of temperature beyond 200°C. This tends to decrease the deformation 
after 200°C and the first-ply failure occurs at 259.5°C, This peculiar variation is due to 
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lig. 3,5 Response of (±45/0/90)2, laminate considering thermal 
dependent ( 1 D) and independent (TID) material properties 


sudden decrease in mechanical properties occurring after the glass transition temperature 
of epoxy at 200* C . from the figure, it is clear that first-ply failure temperature is well 
above the glass transition temperature. Therefore, the response after glass transition 
temperature is insignificant and it can be considered that the attainment of glass 
transition temperature is the failure of the composite. From the literature review, it is 
noticed that there are only few studies which have considered the temperature 
dependency in predicting the responses at elevated temperature. In the subsequent 
analysis, responses with temperature-dependent material properties are restricted to a 
maximum temperature of 200°C. It may be mentioned here that similar response 
characteristics are observed in remaining four laminates listed in Table 3.1. 

3.3.2 Stability envelope 

There could be various combinations of applied temperature and the axial compression at 
which the instability (buckling) may occur. The totality of such critical (bifurcation) 
points in the axial load-temperature space constitutes the stability boundary that separates 
regions of stability and instability. The stability 'boundary defines an implicit relationship 
(or an interaction) between the two parameters; the axial compression and the 
temperature, which, for a (45/-45/0/90)2 S laminate is shown in fig. 3.6. In the interaction 
curve, the axial load N x is nondimensionalised as N x b 2 /E 2 h 3 As can be seen from the 



figure, with an increases in temperature the load required to cause buckling goes on 
decreasing. It is so because the total axial compressive load gets enhanced due to the 
boundary constraints. The response when temperature-independent material properties 
are considered, is very similar to those reported by Noor & Peters (1992) and Argyris & 
Tenek (1995). When thermal dependency is accounted for, the stability boundary moves 
inward as the temperature increases, implying that at a given temperature it requires 
lesser load to cause buckling as compared to the thermally independent case. From the 
figure it is clear that the variation in response (with and without thermal-dependent 
properties) is more pronounced at higher temperatures. 



Temperature, # C 


Fig, 3.6 Stability envelope for ( 1.45/0/90)2, laminate considering thermal dependent (TD) 
and independent (TlD) material properties 

The stability boundaries for all panels listed in Table 3.1 are displayed in one single 
figure 3,7. Although the nature of stability interaction curves is similar in all cases, the 
buckling strength under independent loads (i.e. either compression or temperature) is 
significantly different for different laminates. For example, under the axial load the 
buckling strength is the least for a cross-ply laminate, i.e. (0/90) 4s , and the maximum for 
a (i, 45/0)6 laminate; but the trend changes when only thermal loads are considered. One 
interesting observation is that all interaction curves can, more or less, be approximated 
by straight lines and hence can be represented by a simple equation. 

3.3,2 First-ply failure envelope 

The first-ply failure interaction for a quasi-isotropic laminate, i.e, (45/-45/0/90K under 
combined thermal and inplane compressive loads is shown in fig. 3.8. Such a graph gives 




Fig. 3.7 Stability envelopes considering thermal dependent (TD) and independent (TED) 
material properties 



Fig. 3.8 Stability and failure envelopes for (+45/-45/0/90) 2 , laminate considering thermal 
dependent (TD) and independent (TID) material properties 

information regarding the maximum compressive load that can be applied at a given 
temperature from the consideration of the first-ply failure. The interaction curve shows a 
+ve slope as temperature increases beyond 20°C but it changes at elevated temperatures. 
As is obvious, with temperature rise, the plate can withstand more load than at 20°C 


(zero thermal load). However, at elevated temperatures, the curve is associated with a 
high negative gradient, which implies that the capacity to withstand the axial 
compressive load decreases rapidly. Interestingly, the nature of the first-ply failure 
envelope is distinctly different from that of the stability envelope. In fig. 3.8, the 
inclusion of thermal dependent material properties imparts considerable changes in the 
failure predictions under the combined loading. The graph shows that the plate can 
withstand hardly any loads, beyond 200°C. 



Temperature, 1° C 


(45/-45/0/90) 2i TID 
(45/-45/0/90)2, TD 
(45/-45/0 2 ) 2i TID 
(45/-45/0 2 ) 2j TD 
(45/-45/0 6 ), TID 
(45/-45/0 6 ), TD 
(45/-45) 4s TID 
(45/-45) 4s TD 

(0/90) 4s TID 
(0/90) 4 , TD 


Fig. 3.9 First-ply Failure envelopes considering thermal dependent (TD) 
and independent (TID) material properties 


Figure 3.9 gives failure envelopes for all laminates listed in Table 3.1, with and without 
considering thermal dependent material properties. As can be seen, the peculiar nature of 
the failure interaction curve (as compared to the stability interaction curves) exists for all 
five laminates considered. The interaction curve for the cross-ply laminate, i.e. (0/90) 4s , 
gives a lower bound while that for the quasi-isotropic laminate, i.e. (+45/-45/0/90) 2s , 
yields an upper bound to first-ply failure envelopes. 


In order to gain a better understanding of the peculiar nature of first-ply failure 
envelopes, it is necessary to gain an insight into the deformation response and the 
associated stress-distribution the resultant of which causes the failure. In this regard, a 
detailed study is made for the quasi-isotropic laminate, i.e. ( 45 /- 45 / 0 / 90 ) 2 s, with 
temperature independent material properties. 


Deformation responses in terms of maximum transverse displacement and edge 
displacement are displayed in fig. 3.10 and fig. 3.11, respectively. The transverse 
deflection under temperature rise is already given in fig. 3.5 and it is noticed that 
buckling occurs at 155°C. Figure 3.10 shows the axial load versus the maximum 
deflection variation for a given temperature rise upto the first ply failure. At the room 
temperature (T=20°C) the buckling occurs for a non-dimensionalised value of the axial 
compression equal to 23.5. In addition, the temperature rise is associated with an 
increased transverse deflection and the graph shows that the deflection increases under 
further application of compressive loads for all temperature ranges. 



T = 20° C 
T = 72°C 
T » 122°C 
T-222°C 
T = 371° C. 

T-522 °C 
T - S72 “C 
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Fig.3.10 Transverse deflection for (45/-45/0/90 ) 2j under combined loading at 
different temperature interval 


The inplane displacement is calculated at the center of the edge, x-0 (node 56 in fig. 
3.4). Figures 3.11a & b show the variation of the edge displacement under the axial 
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Fig, 3.11a Edge displacement under axial 
compression at x = 0, node 56 



u (mm) 


Fig. 3.11b Edge displacement under 
temperature rise at x = 0, node 56 


compression and the uniform temperature rise, respectively. From fig. 3.11a it is clear 
that with an increase of compressive load, the edge moves in the +ve x-direction. On the 
other hand, under the temperature rise, the plate exhibits expansion (fig. 3.11b), but 
beyond the buckling temperature, the expansion becomes less prominent followed by an 
inward movement of the edge. Thus, in the pre -buckling range, edge displacements (at 
x~0, a; node=56, 66) under the temperature rise and the axial compression are opposite 
in nature, while in the post-buckling range (near failure) the plate deforms in similar 
ways (edge shortening). It may be noted that edges at y=0, b undergo expansion when 
subjected to the compressive load (in the x-direction) and temperature rise. 

The pattern of the distribution of stresses in the laminate (at the failure location) is of 
paramount importance in the understanding of failure characteristics. For initial 
temperature ranges at which the failure interaction curve shows a higher positive gradient 
(see fig3.8), the failure is noticed in the second layer from the bottom (lamina 
orientation, a = -45° ) at nodes 1 and 121 (diagonally opposite comers). For higher 
temperature ranges, the failure is observed in the 15 th layer (a = -45°; nodes =1, 121). 
The mode of failure in both the case is of transverse type. The stress values are obtained 



at these lamina failure locations. Graphs are plotted between the non-dimensionalised 
axial load and stresses (in GPa) for different initial temperature rise. First, the stress 
distributions at the failure location associated with initial values of temperature under 
combined loading (at node, 1 layer, 2) are cosidered. At this particular location, it is 
observed that the stress cr y has the same values as a x , while other stress components 
(except txy) are close to zero. Figures 3.12a & b give the variation of cj x and while 
fig. 3.12c shows the distribution of stresses o x and txy under the uniform temperature 
rise. Whereas at T=20°C, the plate developes compressive stresses under the axial 
compression (fig. 3.12a) the stress o x is tensile under the temperature rise (fig. 3.12c). 
Thus, at this particular node, the stress o x under two loads is opposite in nature. The same 
holds for the shear stress Txy as well. Thus, the laminate when subjected to a temperaute 
higher than 20°C, is required to be subjected to an increased axial load (compressive) for 
the attainment of first-ply failure conditions. This results in a positive gradient in the 
failure envelope (see fig. 3.8). 



Fig. 3.12a Variation of stress ct x at node 1, layer 2 
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Fig. 3.1 2b Variation of stress t xy at node 1, layer 2 
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Fig, 3.12c Variation of stresses under temperature rise at node 1, layer 2 

Figure 3.13a & b display a x and txy distributions in the 15 th layer at node point 1. At 
higher temperature, failure occurs at this point under combined loads. Here also o x and 
ay have the same values while the stress components x x2 and x yz are neglegibly small. 
Graphs clearly show a substantial variation in stresses before and after buckling. In figs. 
3.13a & b, the pattern of first three curves is significantly different from that of the rest. 




The reason is that the temperature associated with each of these curves is less than the 
critical temperature, whereas for rest of the curves the plate is already in the post- 
buckling state due to the initial temperature rise. It may be noted that the nature of 
stresses at higher axial load is similar to the nature of the stresses under the temperature 
rise (see fig. 3.13a, b & c). It implies that, at the node under reference, at high 
temperatures the plate develops a stress state which is similar to that developed under 
high compressive loads. Because of this reason, the failure interaction curve in fig. 3.8 
shows a steep variation (-ve slope) at higher temperature ranges. 
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Big. 3.13b Variation of stress t*y at node 1, layer 15 
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Fig. 3.13c Variation of stresses under temperature rise at node 1, 
layer 15 

From the above discussion of failure under combined loads, it is clear that it is either the 
axial load (at low value of temperature) or the temperature (at elevated values of 
temperature) which decides the location and the mode of failure. A transition from failure 
due to the axial load to that by the thermal load is occurring at intermediate values of 
temperature and from fig 3.8 it is clear that the failure interaction curve shows a gradual 
transition in the gradient from +ve to -ve. Moreover, the curve shows distinct changes in 
its slope at different intervals. 

As explained earlier, at initial temperatures the laminate failure occurs at comer point 1 
of the 2 nd lamina. During the transition in the gradient in fig. 3.8, the first-ply failure 
occurs at the center of the edges x - 0, a (node 56,66) in laminae 13,15 & 16, the mode 
of failure being shear mode. It is to be noted that a change in slope of the failure 
interaction curve shows a change in failure characteristics and failure locations. It is 
observed that in intermediate temperature ranges, one or more than one of laminae from 
10 to 16 (layer numbers) fails during loading. The distinct slope changes are associated 
with a change in the failure location, lamina and mode of failure. 

3.3.4 Parametric study: effect of number oflayers 

A limited parametric study is carried out by considering the effect of number of layers in 
a symmetric quasi-isotropic laminate, (45/-45/0/90) ns , by varying, n from 1 to 5. The 
analysis is performed by keeping the lamina thickness as constant (thickness of 
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Fig. 3.14 First-ply failure envelopes for (45/-45/0/90) m laminates under combined 
loading. 

individual layer being 0.127 mm). Thermal dependent material properties are taken into 
consideration. First-ply failure interaction curves are obtained and shown in fig. 3.14. 
Graphs are plotted for X* versus temperature (T°C). Here A* is the axial load normalized 
with buckling load for (45/-45/0/90)2s laminate at T=20°C. The maximum temperature of 
200°C is considered as mentioned in section 3.3.1. It is clear that the first-ply failure 
strength of the laminate does not get reduced in the presence of a uniform temperature 
field. 



CHAPTER 4 


CONCLUDING REMARKS 


The present investigation has been concerned with the post-buckling and failure response 
of symmetrically laminated composite plates subjected to combined temperature change 
and in-plane compressive edge load. The simply supported boundary condition, with 
edges free have in-plane movement in the normal direction and restrained against 
tangential displacement, has been adopted same in all cases. Laminates with temperature 
dependent and independent material properties have been analysed. Numerical studies 
have been made of the response characteristics, mainly stability and first-ply failure 
envelopes, of five different types of laminate lay-up that are subjected to combined 
loading. Results have been presented that show the variation in the laminate behavior for 
both temperature dependent and independent cases. An investigation of the stress- 
distributions has been made at laminate failure locations. 

Based on numerical studies, the following observation seem to be justified : 

• Stability interaction curves show a linear variation between the axial compression 
and the temperature rise for all lamination sequences considered. When the 
temperature dependency is accounted for, the laminates require lesser load to cause 
buckling as compared to the thermally independent case. 

• Failure envelopes are of the similar nature for all laminates considered. The 
interaction curves are having a positive slope at smaller values of temperature, while 
at higher temperature ranges the curves are associated with a steep negative variation. 
In both cases, the failure takes place near the comer region, in the transverse mode. 




• In the case of intermediate values of temperatures, the failure envelopes show a 
gradual change in the slope from + ve to -ve. Herein, failure characteristics show 
considerable changes when compared to that at the initial and the higher temperature 
ranges, in regard to the failure location and the nature of failure 

• The failure interaction curve for the cross-ply laminate gives a lower bound, while 
that for the quasi-isotropic laminate yields an upper bound to first-ply failure 
envelopes. 

• When temperature dependent material properties are accounted for, both stability and 
the first-ply failure envelopes (specially the latter) exhibit a distinct change. 

• All the laminates do posses a considerable post-buckling strength in reserve. 
However, since the delamination aspect has not been considered, the decision to use 
this reserve strength would depend upon the maximum permissible transverse 
deflection consideration. 

A few suggestions for further research are : 

• Extensive parametric studies for a better understanding of the response characteristics 
of composite laminates with different boundary conditions and subjected to non- 
uniform temperature fields. 

• Study of effects of geometric discontinuities in the form of cutouts and stiffeners. 

• Enhancing the computational model by incorporating the heat transfer characteristics 
and inelastic material properties at higher temperature. 
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APPENDIX I 


The thermo-mechanical stress-strain relations for the E h lamina of the laminate in the 
principal material directions can be expressed as (Jones, 1975) 
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and as is given in the following section. 


oj ~ C/jSj + C23S2 


where 


c„ - 

IjK J 

1 - V 12 V 21 ’ 

C u 

C 22 " 

e 2 . 

C 33 : 

•-Wl ’ 

C i5 

O t3 ; 

^2 L 
F 

L, 1 

Cj3 

*31 +W32 

_ V 13 + *'l2» / 23 

E 2 E 3 a 

E X E 2 A 

C 23 r ' 

^32 +» / 12 V 31 

_ ^23 +^21 y a 

E t E 3 A 

E X E 2 A 

A = 

l-v n u 2l -V 

23^32 “ V-3, V , 3 - 


_ ^ 12-^2 
1 — V n V 2x 


^12 


C44 ~ G23 


'll y Z2 y \2 


E l E 2 E 3 



and [ Gi, <72, <J3, tj2, T23, tjs f are the stress components and [ej, s 2 , S3, Y12, Y23, Yisf are the 
total strain components referred to the lamina coordinates (1-2-3). The C v are the elastic 
constants of the tf h lamina; E, are the Young's moduli of elasticity in the i th direction, v y 
are the Poisson's ratios that give strain in the f h direction due to stress in the i th direction 
and G 12, G 23 and G 13 are the shear moduli principal plane direction 1-2, 2-3, and 1-3, 
respectively. The cci and thermal expansion coefficients in the material directions 1 
and 2 respectively and AT is the temperature change. 

The plane stress reduced elastic constants in the structural axis (laminate axis), Q ',js of 
the k th lamina, are 


Q11 - Cue 4 + (2Cj2 + ^€33)^ c 2 + C22S 4 

0)2 ~ (S 4 + C 4 ) C 12 + (C)) + C22 ~ 4C3s)^ c 2 

Q13 - (Cu~ C12- 2C33) c 3 s + (C12 - C22 + 2C3s)s 3 c 

Q22 ~ CllS 4 + (2C)2 + 4C33)^C 2 + C22C 4 

Q23 “ (C)i ~ C12- 2C33) s*c + (C12 - C22 + 2C S 3)sc s 

0,66 ~ C)]S 4 + (Cl! - 2C12 + C22 - 2C33. )$ 2 C 2 + (c 4 + 5 4 )C 33 

Q44 = C44C 2 + C55S 2 

Q55 = C44S 2 + C55C 2 


Where O,, 0 Jh c =*cosa, s = sin a and a = lamina orientation which is measured 
positive in the ani-clock wise direction. 



APPENDIX II 


Substituting the eqn. (2.81) in eqn. (2.87) , Kt can be written as 
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It can be shown that 
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for Wj —> nodal degree of freedom, 
for other degrees of freedom 


(A2.2d) 


— L for w, — > nodal degree of freedom. 

(A2.2e) 

0 for other degrees of freedom 
Substituting eqns. (A2.2) in eqn. (A2.1) and rearranging gives the expression for Kt as 
shown in eqn (2.88) 
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